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We present analysis of the impact of structural disorder on the static scattering function of /-armed star 
branched polymers in d dimensions. To this end, we consider the model of a star polymer immersed in a good 
solvent in presence of structural defects, correlated on large distances r according to a power law ~ r~ a . 
In particular, we are interested in the ratio g(f) of scattering intensities of star and linear polymers of the 
same molecular weight, which is a universal experimentally measurable quantity. We apply the direct polymer 
renormalization approach and evaluate results within the double e = A — d, 5 = 4 — a-expansion. We find 
an increase of g(f) with increasing 5. Therefore, an increase of disorder correlations leads to decrease of 
distinction between the size measure of a star and linear polymers of the same molecular weight. 

Key words: polymers, structural disorder, universality, renormalization group 

PACS: 61.41,+e, 61.25.hp, 64.60.ae 



1. Introduction 

Scattering experiments are commonly used in investigations of the structure properties of con- 
densed matter for more than a century (see, e.g. Q). For polymer systems, the quantity of in- 
terest is the static structure function S(k) as function of the wave vector k, representing the 
Fourier transform of the monomer-monomer correlation function 043- The scattering intensity 
I{k) = S(k)/S(0) at small k gives the radius of gyration R g of a single macromolecule: 

J(fc) = l-fc 2 ^ + ..., (1.1) 

where d is the space dimension and (. . .) denotes average over the ensemble of all conformations, 
the macromolecule in a solvent may attain. 

In this paper, we intend to derive some theoretical predictions and quantitative description on 
the peculiarities of scattering experiments with star-like polymers. The star-like polymers are the 
simplest representatives of the class of branched polymer structures, being in a close relationship 
to complex systems such as gel, rubber, micellar and other polymeric and surfactant systems [8|- 
lioj . In particular, some conformational properties of star polymers could be easi ly g eneralized 



to determine the behavior of polymer networks of more complicated structure [111 . Il'2lj . The star 
polymer can be viewed as / linear polymer chains (arms) linked together at the central core (see 
Fig. 1). For / = 1(2) one restores the polymer chain of linear architecture, whereas it has been 
shown that in another limiting situation (/ 3> 1) the star polymer attains the features of soft 
colloidal particle [131 ] . 

A convenient parameter for comparing the size measure of a star consisting of / arms (each of 
lengths N) and a linear polymer chain having the same molecular weight fN is the ratio: 

g(f) ee j^L, (1.2) 

\ gchain/ 
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Figure 1. Schematic presentation of an /-armed star polymer. 



here {R 2 gstBjr ) is the mean square radius of gyration of the star polymer. Note, that the value (|1.2[) 
is an experimentally measurable quantity. Indeed, recalling one has at small k: 

/(fc)star_1 =<?(/) + ..• (1.3) 



^(^)chain 1 

Therefore, comparing the values of the corresponding scattering intensities at the same value of 
the wave vector allows to define the g(f) ratio experimentally. Moreover, it is well established, that 
the gyration radius of a star polymer scales with its total number of segments according to the 
scaling law: 

iK**) ~ V N ) 2v > (1-4) 

with an universal exponent v that also governs the scaling of a single polymer chain of N monomers: 
(^gchain) ~ N 2 " . Therefore, the ratio g(f) of a star and linear chain of the same total molecular 
weight is TV- independent. 

In the pioneering work of Zimm and Stockmayer an estimate for the size ratio g(f) was 
found analytically: 

9(f) = iJ J^- (1-5) 

Inserting / = 1 or / = 2 in this relation, one restores the trivial result g = 1. For any / > 3, ratio 
(|1.5I) is smaller than 1 , reflecting the fact that the size of a branched polymer is always smaller than 
the size of a linear polymer chain of the same molecular weight. Note that in deriving the expression 
(jl.5D the excluded volume effect was neglected (restricting to the idealized Gaussian case, where 
(jl.4l) holds with size exponent v — 1/2). In this limit the result holds for any space dimension. 
Dimensional dependence of the size ratio g(f) is found introducing the concept of excluded volume, 
applied to polymer macromolecules by P. Flory. It refers to the idea that any segment (monomer) 
of macromolecule may not occupy space that is simultaneously occupied by another segment; this 
causes a swelling of a polymer in solution with size exponent v{d) > 1/2. Later analytical |15l4l7| 
and numerical [18l422| studies have found the value of g(f) to increase if the excluded volume effect 
is taken into consideration. 

While reliable estimates are known for the value of g(f) for polymers in a good solvent, there 
are no similar estimates for the case when the polymer is immersed in a good solvent at presence 
of structural impurities or a porous environment. However, such estimates are of a great impor- 
tance for an understanding of the behavior of macromolecules in colloidal solutions [23[ or near 



microporous membranes [24[. The density fluctuations of such obstacles lead to a large spatial 
inhomogeneities, which often produce pore-like fractal structures (25[. Such a disordered (porous) 
environment may be found, in particular, in a biological cell, composed of many different kinds of 
biochemical species (26[. It has been proven both analytically 27] and numerically 28], that the 
presence of uncorrelated point-like defects of weak concentration does not change the universality 
class of polymers. Here, however, we address a case, where the structural obstacles of the envi- 
ronment are spatially correlated on a mesoscopic scale j29[. Following Ref. [3(j, this case may be 
described by assuming the defects to be correlated on large distances r according to a power law 
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resulting in a pair correlation function: 

h(r)~ r - a . (1.6) 

For a < d, specific situations giving rise to such a correlation function include defects extended in 
space e.g. the cases a = d — 1 (a = d — 2) may be represented by lines (planes) of defects of random 
orientation, whereas non- integer values of a may include obstacles of fractal structures (see [3(J EI 
for further details). The impact of long-range-correlated disorder on the scaling of linear and star 
branched polymers has been analyzed in previous works (3l| by means of the field-theoretical 
renormalization group (RG) approach. The aim of the present paper is to evaluate analytically an 
experimentally measurable ratio (|1.3[) for a star polymer in a solvent in the presence of structural 
defects correlated according to (|1.6p . 

The layout of the reminder of the paper is as follows. In the next section, we develop a description 
of the problem in the frames of the Edwards continuous chain model. In section III, the direct 
polymer renormalization method is briefly described; the results of its application to the problem 
under consideration are presented in Section IV. Conclusions and an outlook are given in Section 
V. 



2. The model 



We consider star polymers with / arms in a solution in the presence of structural obstacles. In 
the frames of the Edwards continuous chain model [l| , each arm of the star is presented by a path 
rj(s), parameterized by < s < Si, i = 1, 2, . . . , /. The central branching point of the star is fixed 
in space, so that: r*i(0) = ... = r}(0) = 0. We take the contour length of all arms to be equal: 
S± = . . . = Sf = S. The partition function of the system reads: 



Z f (S)= / 2>{r}Ip(*(0))«p 



f ° 



dri(s) 
ds 



f f 

| E/ ds 'J ds"6(r t (s")-r 3 (s>))+J2 J dsV(f t (s)) 



(2.1) 



Here, a multiple path integral is performed for the paths r*i, . . . , rf and the product of (^-functions 
reflects the star-like configuration of / chains. The first term in the exponent represents the chain 
connectivity, the second term describes the short range excluded volume interaction with bare 
coupling constant bo, and the last term contains random potential V(fi(s)) arising due to the 
presence of structural disorder. Denoting by (. . .) the average over different realizations of disorder, 
the first moment of the distribution is: 



V(f(s)) = p 

with po being the density of obstacles. Let us introduce the notation for the second moment: 

Viris^Viris')) = h(r(s) - f(s')). (2.2) 

Note that dealing with systems that display randomness of structure, one usually encounters 
two types of ensemble averaging, treated as quenched and annealed disorder (32|. The annealed 
case amounts to averaging the partition sum of the system over the random variables, whereas in 
the quenched case the free energy (or the logarithm of the partition sum) is to be averaged; the 
replica formalism is usually exploited in the last situation. In principle, the behavior of systems 
with quenched and annealed disorder is quite different. However, as has been shown in a number 
of works (33l - [3^ |. the distinction between quenched and annealed averages for an infinitely long 
single polymer chain is negligible, and in performing analytical calculations for quenched polymer 
systems one may thus restrict the problem to the simpler case of annealed averaging. To average 
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the partition function of the system over different realizations of obstacles, we make use of the 
relation: 



exp 



where M n (x) are nth cumulants of the random variable x: Mi(x) — x, M%{x) — (x — x) 2 etc. 
Noticing that only the last term in (|2.1I) contains random variables and taking into account (|2.2p , 
we obtain: 

r I 

Z f (S)= DWl]^(0))e" fl - (2.4) 
J 1=1 

with an effective Hamiltonian: 

y S 2 f & S 



1=1 = 1 

f s s 

\il J ds 'J ^s"h(n{ S ")^r 3 { S '))-p a fS- l -plfS 2 . (2.5) 



2 

The last two terms in f|2.5[) correspond to a trivial constant shift which will be omitted in the 
following analysis. Note also, that in (|2.5[) we do not take into account terms generated by higher- 
order correlations of the type (I2.2p . for the problem under consideration these are irrelevant in the 
renormalization group sense. 

The case of structural disorder in the form of point-like uncorrelated defects corresponds to 
h(r{s") — r(s')) = vqS^s 1 ') — r(s')) where vq is some constant. One immediately reveals that in 
this case one can adsorb the effect of disorder into the excluded volume coupling constant passing 
to the coupling: bo = bo — vq. This conclusion was obtained for the case of polymers in quenched 
disorder by Kim [27| on the basis of a refined field-theoretical study; in the present case of annealed 
disorder this is a straightforward result. 

We address the model where the structural obstacles are spatially correlated at large distances 
r according to (|f .6p . Taking into account that the Fourier transform of the correlation function at 
small k is related to its large-r behaviour via: 

= W) -f j (s')\- a ^w [ Akk a - d e k ^ s "^^\ (2.6) 



one is left with a model with two couplings bo and wo- Note that coupling bo must be positive, 
which corresponds to an effective mutual repulsion of the monomers due to the excluded volume 
effect. The coupling wq is positive as results from the Fourier image of the correlation function. 

Performing dimensional analysis for the terms in (|2.5[) one finds the dimensions of the couplings 
in terms of dimension of contour length S: [bo] = [S] db ° , [wo] — [S] d ™° with d(, = (4 — d)/2, 
d W() = (4 — a)/2. The "upper critical" values of the space dimension (d c — 4) and the correlation 
parameter (a c — 4), at which the couplings are dimensionless, play an important role in the 
renormalization scheme, as outlined below. 



3. The method 

To study the universal properties of polymer macromolecules in solutions, it is convenient to 
apply the direct renormalization method, as developed by des Cloizeaux The efficiency of this 
approach comes on the one hand from its close relation to the concepts of field theory [38j, and 
on another hand from providing a considerably simpler treatment of a variety of complex polymer 
systems. 
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In the asymptotic limit of an infinite linear measure of the continuous polymer curve (cor- 
responding to an infinite number of configurations), one observes various divergences. All these 
divergences can be eliminated by introducing corresponding renormalization factors, directly asso- 
ciated with physical quantities. This postulates the existence of a limiting theory, describing sets 
of very long polymers. 

As a first step within this theory, we consider the size measure of an /-arm star polymer, 
given by the mean square end-to-end distance of its individual arm. When evaluated in terms of a 
perturbation theory series in bare coupling constants {Ao}, this reads: 

(Rl) = ((r(S) - f(0)) 2 } = Xo({A })S. (3.1) 

Here, the averaging is performed with respect to a corresponding effective Hamiltonian, and 
Xo({An}) is the so-called swelling factor, reflecting the impact of interactions on the effective size of 
the macromolecules. For the case of a Gaussian chain (all couplings Ao = 0) one has Xo({0}) = 1. 
Recalling the scaling of the polymer size with its molecular weight: 

(R 2 e ) ~ N 2 " ~ S 2v , (3.2) 

one finds an estimate for the effective critical exponent ^({Ao}): 

2^({A })-f-5— ln Xo ({A }. (3.3) 



The second renormalization factor Xi({Ao}) is introduced via: 

= Lxi({Ao})] 2 . (3.4) 



Zf(S) _ r .. n,2 
Z%S) 



Here, Zf(S) is the partition function of an /-arm star polymer and Z®(S) - the partition function 
of an idealized Gaussian model. It is established, that the number of all possible conformations of 
an /-armed star polymer scales with the weight of the macromolecule parametrised by S as: 

Zf{S) ~ ^ s {fSyt-\ (3.5) 

Here, the 7/ are additional universal critical exponents, depending only on the space dimension d 
and the number of arms / (exponents 71 = 72 = 7 restore the value for the single polymer chain) , 
[i is a non-universal fugacity. In similar way as for the size measure, one finds from the scaling 
assumption (|3.5p an estimate for an effective critical exponent 7/({Ao}), governing the scaling 
behavior of the number of possible configurations as: 

7/({Ao})-I _ ainxi({A }) 

2 ~ 5 ds ' [6 - b) 

The critical exponents (|3.3[) and (|3.6p . presented in the form of series expansions in the coupling 
constants {Ao} are, however, divergent in the asymptotic limit of large 5*. To eliminate these 
divergences, the renormalization of the coupling constants is performed. Subsequently, the critical 
exponents attain finite values when evaluated at the stable fixed point (FP) of the renormalization 
group transformation. Note that the FP coordinates are universal, in particular the scaling of a 
single polymer chain and that of a polymer star is governed by the same unique FP. Therefore, to 
evaluate the FP coordinates in the following analysis we restrict ourselves to the simpler case of a 
single chain polymers (/ — I). To define the coupling constant renormalization, one considers the 
second virial coefficient of a polymer solution, given by the relation: 

here, n is the osmotic pressure, /? = 1/ksT, C is the number of monomers per unit volume, 
and Z\(S) is the partition function of a single polymer chain. Z\ (S,S) are contributions into 
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partition function of two interacting chain polymers, having dimensions Z\ (S,S) ~ [S'] 2 [Ao]. The 
renormalized coupling constants A^ are thus defined by: 

Afl({A }) - -[x 1 ({A })]- 4 Z Ao (L,L)[2 7 r Xo ({Ao}i]- (2 - d ^), (3.8) 

therefore: 

n/ 3 = C + \ E ^RC 2 [2n X o({Xo}Lf- d ^ + . . . , (3.9) 



In the limit of infinite linear size of the macromolecules the renormalized theory remains finite, 
such that: 

lim A H ({A }) = X%. (3.10) 

Moreover, for negative d\ < the macromolecules are expected to behave like Gaussian chains 
in spite of the interactions between monomers, thus each AJj = for corresponding d\ < 0. It 
is therefore proper to choose {\r} as expansion parameters, which remain finite for S — >• oo and 
which are also rather small close to the critical dimensions of the corresponding couplings. The 
concept of expansion around in small deviations from the upper critical dimensions of the coupling 
constants, thus naturally arises. 

The flows of the renormalized coupling constants are governed by functions f3\ R : 

,0A fl ({Ao}) 



dS 



(3-11) 



Reexpressing {Ao} in terms of renormalized couplings A# according to (|3.8[) . consequently the 
fixed points of the renormalization group transformations are given by the common zeros of the 
/3-functions. Stable fixed points govern the asymptotical scaling properties of macromolecules in 
solutions and allow e.g. to obtain reliable asymptotical values of the critical exponents p. 31) and 
(EH). 

4. Results 



We start by evaluating the partition function (|2.4p of the model with an effective Hamiltonian 
(|2.5I) . performing an expansion in coupling constants bo, wq: 



Zf{S) 



Vre i=1 ° v 



1 - T E / ds 7 ds "l dfce* fc >^">-^ s ')) 



+ T E J ds ' j ds " J dkk^e^^-^')) +...). (4.1 ) 

i,i=l o 

here the Fourier-transform of 5-function is exploited and the last term originates from the Fourier 
transform of the function h at small k (see Eq. (|2.6[1 ). Below, we will consider the one- loop approx- 
imation, keeping only the first-order terms in bo, wq in the expansions. One may rewrite: 



(4.2) 



ds 



(4.3) 



taking into account that rj(0) = rj(0) — in our model of a star-shaped polymer. Making use of 
the identity: 



exp 



_i /WdnM 
2./ ds 



+ ik 



dfj(s) 
ds 



ds 



exp 



drj(s) 
ds 



ik 



(4.4) 
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and taking into account that: 



■DC 

I 



dxe 



(4.5) 



we receive: 



Z f (S) = Z%S) 



1-6o(2tt)-3 [ / / ds" J ds'(s" -s')- d/2 + 





s s 

/(/ ~ 1} J ds" J ds'(s' + S ")- d ' 2 \ + 



s 

<r,A2:>y- { f J As" J d.s'(.s" - .«'r" " t 





I (Jl „'\— a/2 i /(/ 1) 



s s 
ds" fds'(s' + s")- a / 2 







.(4.6) 



In the last equation, the Gaussian integration over k is performed and the notation Z®(S) is 
introduced for the partition function of "unperturbed" Gaussian model: 



a (=10 V ' 



Z%S) = J Pre a *=i6 v ds 7 . (4.7) 

In what follows, we will use the diagrammatic representation of the perturbation theory series (see 
Fig. 2). Performing the integrals in (|4.6p and introducing dimensionless couplings 



6 = 6 (2^)- d / 2 S 2 - d / 2 , w = w (2ir)~ a / 2 S 2 - a / 2 



(4.8) 



we obtain: 



Z f (S) = Z°(5) 
4w 



4/; 



(2-d)(4-d) 
( /+ /C/^ (2 W»_ a) ) 



(f+ f ( f 2 1) (2 2 - rf / 2 -2)j 



(2-a)(4- 

Finally, one may perform a double e = 4 — <i, <5 = 4 — a-expansions: 

7(3-/) , /(3-/) , /(/-l) 



(4.9) 



2/(5) = 1 + 6, , 

/ /(3-J) /(3-/) /(/-l) 
\ J 2 2 



ln(2) 



ln(2) 



(4.10) 




Figure 2. Diagram contributions to the partition function of a 4-arm star up to the first order in 
the coupling constants. Dotted line denotes possible interactions between points s' , s" , governed 
by couplings bo and wo- Integrations are to be performed over all positions of the segment end 
points i.e. over all mutual interaction points within a single arm and between different arms. 



xxxxx-7 



V. Blavatska, C. von Ferber, Yu. Holovatch 



The averaged squared end-to-end distance (i? 2 ) of a single arm of a star polymer may be 
calculated using the identity: 

_d 

' dq 

where 



(Rl) = ((f(S) - r(0)) 2 ) = -2d^(e^(r(s)- m)) \ q=Q (4 . n) 



fVre-^(---)ULmO)) 



<••■>-' — wr ; • ( " 2) 

Following the same scheme as described above for the partition function, we find: 

46 4w 



( ^ = Sd { 1+ (4-d)(6-d) ~ (4-a)(6-,) J ' (4 ' 13) 
We may therefore define a swelling factor Xo{bo,wo) (cf. (|3.ip ) as: 

»("■-» " (' + (4-^6-0 - (4-.K6-.) ) = (' + J< 2 " " " f < 2 " ") ' < 4 ' 14 ' 

Now we return to the calculation of the gyration radius. The gyration radius of a star polymer 
in a solvent in the presence of correlated defects may be written as: 



/ / 



(/star 





We rewrite: 



^ £ (^(-2) -fy(«i)) 2 ^ = -2d^(e^-(-W-^(-)))| g=0 , (4.16) 

where d is the space dimension. First let us consider the zero-loop order of the expansion of (|4.16[) 
in coupling constants. A diagrammatic representation is given in Fig. 3. The analytic expression, 
corresponding to the diagram (a) reads: 



f(e i $(?iM-M*t))\ = /e~^" (s2 ~ Sl) , (4-17) 

whereas the diagram (b) gives: 

Note, that here and in what follows the Gaussian integral in (|4.7p cancels against the numerator, 
when averaging (|4.12l) is performed. One thus distinguishes between two types of contributions: 




Figure 3. Diagrammatic presentation of contributions into radius of gyration of star polymer in 
zeros order of perturbation theory. 
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Figure 4. Diagram contribution into the gyration radius at the one-loop level. 



one resulting from insertions si, s 2 along the same arm of the star, and the second corresponding to 
insertions located on two different arms. Taking the derivative with respect to g, and evaluating for 
q = according to (|4.16l) , we have for the radius of gyration of a star polymer in the unperturbed 
(Gaussian) case: 



C^gstar) 



(fsy 



f I ds 2 / dsi (s 2 - s x ) 

o 



/(/-I) 



dsi / ds 2 (si + s 2 ) 



dS 3/ - 2 



(4.19) 

Now, let us perform calculations up to the first order of the perturbation theory expansion 
in the couplings bo, wq. In Fig. 4, we present diagram contributions to (|4.16D . see Appendix for 
details. We have: 



{j^ (^(s 2 )-r;( S i)) 2 ^-2d^[/(/ 1 (d)+/ 1 (a)+/ 2 (d)+/ 2 (a)+/ 3 (d) + 

J 3 (o) + h(d) + h(a) + h(d) + h(a) + 7 6 (d) + J 6 (o) + h{d) + 7 7 (a)) + 
f(f-l) 

2 (W) + h(a) + h{d) + 7 9 (o) + ho(d) + 7 10 (a) + I u (d) + 7 n (a) + 
/12(d) + 7 12 (a) + /13(d) + /13(a) + /14(d) + I 14 (a) + I 15 (d) + /15(a) + 
7 16 (d) + 7 16 (a) + I 17 (d) + 7 17 (a)) + /(/ ~ ^ (/ ~ 2) (/ 18 (d) + 7 18 (a) + 

7 19 (d) + /i 9 (a)/ 20 (d) + 7 20 (a)) + /(/ ~ 1)(/ 24 2)(/ ~ 3) (/ 21 (d) + 7 21 (a))] | g = . (4.20) 

Here, /i(d), 7(a) are integrals listed in the Appendix. Performing the integration according to 
(|4.15p and evaluating the double e, 5-expansion we find for the gyration radius: 



dS 
6? 



(^ S tar) = ^(3/-2) 



1 



(b-w) 



13 13 (/-!)(/ -2) 4(/-l)(3/-5) 
12 2 3/ - 2 nlj 3/ -2 



(4.21) 
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The result for the radius of gyration of a single chain of total length fS is straightforward: 



(R 2 



lin) 



d(Sf) 



pchain/ g 

d(Sf) 



l + bf 



e/2 



13 
12 



6 



13 



12 



- ™f /2 



13 



13 
12 



1 + )-v> T-T7T )+(b-w)Hf) 



12 



(4.22) 



Finally, we need to perform the renormalization of coupling constants. To this end, we need 
the contributions to the partition function Z(S, S) of the system of two interacting polymer chains 
of the same length L. In the diagrammatic representation of this function, one takes only those 
terms into account, which contain at least one interaction line (see Fig. 5). In general, performing 
a thorough dimensional analysis of the contributions, produced by different diagrams, we find two 
distinct classes of the diagrams. The first class of graphs produces terms which behave as [S 1 ] 2 ^], 
the sum of all such terms contributes to the function denoted by Zi, (S, S) ; the second class of 
diagrams behaves as [S] 2 [u>o] and thus contributes to the function Z Wo (S,S). 



We find the contributions to functions Zb g (S, S) and Z WQ (S, S) as: 



Z bo (S, S) 



326 



-b S 2 

2i-d/2 



1 



8w 



(2-d)(4-d) (2-a)(4-a) 



10 



(d-2)(d-4)(d-6)(d-8) 



32 



w 



2a-d + 2 4 -< 2a - d )/ 2 -10 



64u> 



b (2a -d- 2) (2a 

a + 2 4 -"/ 2 - 



d-4)(2a 
10 



(a-2)(a-4)(a-6)(a-8) 
Z Wo (S, S) = w S 2 



d-6)(2a-d-8) 



8w 



(2-d)(4-d) (2-a)(4-a) 
We may therefore define the dimensionless renormalized coupling constants wr by (|3.8[) 



(4.23) 



b R = 



Wr 



(Z(S)) 



■ Z " g 'f 2 ^ [2-KXo(bo,wo)S]~ 



(4.24) 



(Z(S)) 

The RG flows of the renormalized coupling constants are governed by corresponding /3-functions: 



P bR (b,w)=2S 



dbR 
dS 



/3 WR (b,w) = 2S 



dw R 
dS ■ 



(4.25) 



(a) 



(b) 



(c) 



(d) 



Figure 5. Diagrammatic presentation of contributions into the partition function of two inter- 
acting chains. 
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Reexpressing in (I4.25P b and w in terms of the renormalized couplings (|4.24|) . we finally have: 



0b R (b R ,w R ) = (4 - d)b R - b\- — - + w R b R - + 326 



> 2d , 2d „ d + 2 4 - d / 2 -10 

^Trf + w * b n— a + 32 ^ (d _ 2)(d _ 6)(d _; 



2 2a-d + 2 4 -( 2a - d )/ 2 -10 , a + 2 4 - a / 2 -10 



(2a-d-2)(2a-d-6)(2a-d-8) (o - 2)(a - 6)(a- 8)' 

Pw R (b R ,w R ) = -(4 - a)wfl - w 2 R -^— + b R w R -^— (4.26) 

d — a o — a 

Performing a double e, 5 expansion and keeping terms up to linear in these parameters, we then 
have the RG functions: 

Pb = eb R - 8b R - Aw R + 12b R w R , 

Pw = -Sw R - 4w 2 R + 4b R w R . (4.27) 

The fixed points b* R , w R of the renormalization group transformations are defined as the common 
zeros of the RG functions (|4.27l) . We find three distinct fixed points that determine the scaling 
behavior of the system in various regions of the a and d plane: 

Gaussian : b* R = 0, w* R = 0, stable at e, 5 < (4.28) 

Pure: b* R = ^, w* R = 0, stable at 5<e/2, (4.29) 

8 

LR: b* R = lf S n ,w* R = 5 ^~ 26 ) , stable at e/2 < S < e. (4.30) 

Here and below the index LR means that corresponding quantity is evaluated in the region of d, a 
plane, where the influence of long-range-correlated disorder is relevant. 

We can also find estimates for critical exponents, governing the scaling of star polymers in solu- 
tions in the presence of correlated disorder. Making use of definition (|3.3[) . recalling the expression 
or renormalized scale (|4.14p and passing to the renormalized couplings we have: 

v{b Rl w R ) = - + - — . (4.31) 

Evaluating this expression at fixed points of the renormalization group transformation (|4.28l) - 
(|4.30p . we find the corresponding estimates for the size exponents: 

^Gaussian _ ^ ^ 

^ mc = bh (4 - 33) 

^ = \ + { ( 4 -34) 

Note that these exponents govern the scaling behavior of macromolecules in the regions of stability 
of the corresponding fixed points (|4.28[) - (|4.30[) . 

Similarly, evaluating (|3.6p taking into account (|4.10l) one finds in terms of renormalized cou- 
plings: 

„ x 1 , b R f(3-f) w R f(3 - f) 
l(b R ,w R ) = l + . (4.35) 

Again, evaluating this expression at the fixed points (|4.28p - (|4.30p . we find: 

Russian = ^ (4.36) 

7/ Pure = ^ + ^/(3-/), (4.37) 
7/ R =5 + 5/(3-/)- (4-38) 
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At / = 1 one restores corresponding exponents for a single polymer chain. Note, that the first-order 
expressions for the fixed point coordinates (|4.28[) - (I4.30[) and for the critical exponents (|4.32j) - (j4.34j) 
and (|4.36[) - (|4.38[) were obtained here for the annealed system. Thus, with ^ LR and 7^ R as obtained 
in the regime of annealed disorder, we restore the corresponding exponents, governing the scaling 
of polymers in solutions in presence of quenched long-ranged correlated disorder, studied by us 
previously [31[. This supports the statement of equivalence between the annealed and quenched 
averaging when dealing with polymer systems. 

Finally, we obtain estimates for the size ratio g 
Km and KM- 



of star and linear polymers, recalling 



9(f) 



7 Gauss 
y ) 
^Purc 



,LR 



e, 8 < 0, 
6 < e/2, 
e/2<5< e. 



with: 



Gauss 



3/ -2 



P 

!■•„,, _ 3/ — 2 
P 

LK = 3/-2 

p 



1 - 

1 - 



e /13 (/-!)(/ -2) 
8 V 2 3/ - 2 
5 (13 (/-!)(/ -2) 
4 V 2 3/ - 2 



,„ (2) 4(/-l)(3/-5) +ln(/) 

-M2) 4(/ -^- 5) + M/) 



(4.39) 

(4.40) 
, (4.41) 
(4.42) 



With g Puie we restore the size ratio for the case of polymers in pure solvent, evaluated previously 
in Refs. El, El. 



5. Conclusions and Outlook 



In the present paper, we have studied the impact of structural disorder on the static scattering 
function of /-arm star branched polymers in d dimensions. To this end, we consider the model 
of a star polymer immersed in a good solvent in presence of structural defects, correlated at 
large distances according to (ll.6[) with parameter a [30( . The impact of such long-range-correlated 
disorder on the scaling of linear and star branched polymers has been analyzed in our previous 
works [311 ]. Here, we are interested, in particular, in the ratio g(f) of scattering intensities of star 
and linear polymers of the same molecular weight, which is a universal experimentally measurable 
quantity (see (jl.3p ). We applied the direct polymer renormalization approach [l[ and evaluate 
results within the double e = 4 — d, 5 = A — a-expansion. 

Let us analyze the expressions obtained for the size ratio of star and linear polymers in solution 
in the presence of structural defects. First of all, as far as S = 4 — a-prefactor in (|4.42[) is positive for 
all / > 2, one immediately concludes, that the stronger the correlations of disorder (i.e. the larger 
is parameter 6), the larger is the size ratio g(f) and thus, the smaller is the distinction between 
the size measure of a star and linear polymers of the same molecular weight. In Fig.[S^, we present 
estimates for the size ratio (|4.42D as a function of number of arms /, obtained by direct evaluation 
at several fixed values of the correlation parameter 5. Besides an expected increase of g LR (/) with 
growing 5 at each fixed /, we notice also a decrease of the size ratio with / for any fixed S. This 
can correspond to a shrinking of the effective radius of star polymers with increasing functionality. 

Another interesting aspect results from comparing qualitatively the impact of long-range- 
correlated disorder on the size ratio with the impact of the excluded volume effect in pure solvent. 
We consider the ratio: 



g Purc (/) 

3 LR (/) 



l 



13 (/-!)(/ -2) 
2 3/ - 2 



ln(2) 



4(/-l)(3/-5) 
3/ -2 



Hf) 



(5.1) 



In Fig. we present the evaluation of this ratio as function of / at three dimensions (taking 
e = 1) and several fixed values of correlation parameter S. Note, that correlated disorder with 
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/ 
(a) 



1.4 




0.6 ' ' ' ' ' ' ' 1 

23456789 10 

/ 
(b) 



Figure 6. Left: the size ratio (|4.42|l as a function of / at different values of correlation parameter 
8. From below: 8 = 0.1, 8 — 0.5, 8 — 1.0. Right: the ratio (|5.ip as a function of / at different 
values of correlation parameter 8. From above: 8 = 0.1, 8 — 0.5, 8 = 1.0. In both (a) and (b) we 
fix the value of parameter e — 1. 



5 = 0.5 plays the role of "marginal", separating a region where the ratio (|5.1[) increases with / (at 
5 < 0.5) and region where it decreases with / (any S > 0.5). 

Let us recall, that with the critical exponents f LR (|4.34[) and 7y R (|4.38l) that we obtained in 
the present paper in the regime of annealed disorder, we restore the corresponding exponents, gov- 
erning the scaling of polymers in solutions in presence of quenched long-ranged correlated disorder, 
studied by us previously [3lj . This supports the statement of equivalence between the annealed 
and quenched averaging when dealing with polymer systems. Thus our qualitative estimates for 
the size ratio of star and linear polymers in annealed correlated disorder (|4.42l) also holds for the 
case of quenched systems. 

Note that our results are based on the first-order perturbation theory expansions and give 
rather qualitative description of the impact of structural disorder on the quantities of interest. To 
obtain the reliable quanitative estimates, the higher order analysis would be worthwhile, which is 
the subject of our forthcoming studies. 
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Appendix 

Here, as an example we evaluate the analytic expression corresponding to diagram (9) in Fig. 4 
(shown in more details in Fig. 7), presenting contributions into radius of gyration of star polymer. 
We have for this contribution: 



( e -«(ft(« a )-ft(«0)) (9) =-b / Vre a *=i6 V ds ' e-^M-^W) x 



(9) 

s s 

dfee - ^^ 8 ')--'^^'')) 



S 2 o 
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Rewriting the last exponent: 

- ik(h(s') - f 2 (s"j) = -ik((n(s') - fi(s 2 )) + (n(s 2 ) - r!( Sl )) + (fi(si) - r 2 (s"))), (5.2) 
and making use of (|4.3D - (|4.5[) . one arrives at: 

/ s 

h(d) 



_I £ fds( dr * M " 

( e -Wi(s2)-ri( Sl ))\ _ / x)re <=1 ° v ds ' 



(5.3) 



with: 



s s 

J 9 (d) = J ds" J ds' J dke-^ <- s ■■)- AS ^(^- S i) = 

s 2 

s s 

^ (2n)- d ^e-^^- s ^ J ds" J ds'(s" + s')- d/2 e^ L ^T 

S 2 

Taking the derivative over q according to (|4.16p we receive: 

S S 

^^\ q ^ = (2n)-^d(s 2 - Sl ) J ds" J ds'(s" + sT d/2 - 

s 2 

S S 

-(2TT)- d / 2 d( S2 - Sl ) 2 J ds" J ds' {s" + S')- 1 ' 2 - 1 = 
S 2 

/2 d{s 2 - si) 



(5.4) 



(27r)- d / 

-(27T)- rf / 2 



(l-d/2)(2-d/2) 
rf( S2 -si) 2 



(25) l-d/2 _ {S + S2 )l-rf/2 _ 5 l-d/2 + s l-d/2l _ (5 5) 



(1 - d/2)(-d/2) 

Finally, the contribution into the gyration radius (|4.15[) can be found by integrating over si, s 2 . 
Note, that the Gaussian integral in (I5.3[) cancels against the numerator, when averaging (|4.12l) is 
performed. 

The expressions corresponding to other diagrams in Fig. 4 are listed below (note, that 1(d) and 
1(a) arise, correspondingly, when treating interactions with couplings bo and wq respectively). The 
factors (2n)~ d / 2 and (27r)~ a / 2 in front of each integral are omitted . 



S S2 



„2 (. a -,')2 



h(d) = c-V(^i) / ds" / ds' (s" - s ')- d / 2 e-T^T, 



S2 SI 




Figure 7. Example of diagrammatic contribution into the radius of gyration of star polymer 
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S2 Si 

2 l2 



h(d) =e-^( S2 " Sl ) y da'^ daV'-s 7 ) - ^ 2 ^!^ 

si 
S si 

J 3 (d) = e -^^- Sl ) y ds'Y da' (a" - a') 



I 1 (°2-»l) Z 
d / 2 e T (o"-a') 



S 2 



J 4 (d) = e"£ (s2 ~ Sl) y da" y da' (a" - a')" 



d/2 



S2 S2 



Si S 

-d/2 



<i/2 e V( s "- s i) 



J B (d) =e-^( S2 - Sl ) y da'^ da' (a" - a') 
o o 

// 

/ 6 (d) = e - £ (S2-Sl) y da" y da' (a" - a') 

si si 

S si 

/ 7 (d) = e-^ S2 " Sl ) Jds" Jds' (a" + a')- d : 
o o 

S S 2 

I s (d) = e-£^- s ^ J As" J da'(a" + a')- d/2 e^WT 

si 

s s 



,,„ ^ ('2-»ir 

d / 2 Q~ {s' + s'l) 



h{d) = c -^(^-s!) y ds "y ds / (s" +s ') 

s 2 

s s 

7 10 (d)=c-^( S2 - Sl ) y da'^ da'(a" + a')- d/2 , 
o o 

S s" 

/n(d) = e -^( S2 - Sl ) y da'^ da'(a"-a')- d/2 , 
o o 

s 

/12(d) = C -^( S2 +^) y da'^ da' (a" - a')- d / 2 , 

Sl Si 

S 81 ,2 

/13(d) =c-^ S2+Sl ) y da'^ da'(a" -a')- d / 2 c^T^T, 

si 

Si 8 

/14(d) = e -^( S2+Sl ) y da'^ da' (a"-a')- d / 2 c^( s "- s '), 

s s 

J 1B (d) = e~£ (S2+Sl) y da'^ da'(a" + a')- d / 2 e^T^T?T, 



S2 Si 
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S2 S ^ 

he{d) = c -^+ s i) f ds" f ds'{s" + s')-^ 2 e^ i ^fvr. 



si 

S2 



j 17 (d) = e -£ y d^y d S '( S "+ S ')- d/2 c^( s " +s ' 



S s 2 

z 18 (d) =e -4(«*+«o y d^yds' («"+*' 



Cf2)j 





s s 

/ 19 (d) =e -^( S2+Sl ) y ds'^ ds' {s" + s')- d ' 2 e^T^-) . 

S 2 o 

= c -^ s2+si > y d^y d S ' (s" - s')-^ 2 , 


s s 

hi(d) = e-^ {82+Sl) [ds" Jds' (s" + s')- d/2 - 
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Biuihb Gesjia^y Ha cTaTiiHHy cpyHKniio po3ciaHHH 3ipKOBiix nojiiMepiB 
B. BjiaBaiibKiP, K. cpOH Oep6erj2^1, K). TojiOBa^P 

IHCTHTyT <pl3HKH KOHfl6HCOB£LHHX CHCT6M HAH y KpaiHH , 79011 JlBBiB, YKpaiHa 

/locjiiflHHij,BKHH ueHTp npiiKjia,z],HOi MaTeiviaTHKH, ymBepciiTeT KoBeHTpi, CV1 5FB KoBeHTpi, 
BejiiiKo6pHTaHm 

IHCTHTyT TeopeTHHHo'i <pi3HKH II, ymBepcHTeT im. reHpixa TeiiHe, D-40225 /lioccejiBflopcp, HiMeHHima 

npe^CTaBjieHO aHajii3 BnjiHBy CTpyKTypHoro 6e3jia,zry Ha CTaTHHHy cpyHKuiio po3ciHHHH /-rijiKOBoro 
3ipKOBoro nojriMepa y rf-BHMipHOMy npocTopi. Po3rji5myTO Mo^ejiB 3ipKOBoro nojiiMepa y xoponiOMy 

p03MHHHHKy B npHCyTHOCTi CTpyKTypHHX fleCpeKTiB, CKOpejlBOBaHHX Ha B6J1HKHX Bi^ajIHX r SriflHO 

CTeneHeBoro 3aKOHy ~ r~ a . 3oKpeMa, mh HiKaBHiviocB BiflHOineHHHM g(f) iHTeHCHBHOCTen po3cmHH5i 
3ipKOBoro Ta jiimiiHoro nojiiMepiB OflHaKOBoi MOjieKyjinpHOi Macn, mo e ymBepcajiBHOio, eKcnepHMeHTajiBHO 
cnocTepeacyBaHOio bg jihhhhoio . Mh 3acTOCOByeMO mctoa npaMoro nojiiiviepHoro nepeHopMyBaHHH i 
BHKopHCTOByeMO noflBiiiHHH e = 4 — d, 5 = 4 — a-po3Kjia,z],. 3HanfleHO 3pocTaHmi bcjihhhhh g(f) 13 
3pocTaHHHM napaMeTpa S. TaKHM hhhom, 3pocTaHHH KopejiHuiii 6e3jia,n,y npHBOflHTb flo 3MeHHieHHH 

Bi^MiHHOCTi Mi^K p03MipOM 3ipKOBHX Ta JlimHHHX nojiiMepiB i3 OflHaKOBOK) MOJieKVJIHpHOK) Baroio. 

KjnoHOBi cjiOBa: nojiiiwepH, CTpyKTypHHH 6e3jiafl, yHiBepcajiBHicTB, peHopMajii3au,iitHa rpyna 
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